BCS-BEC crossover in a strongly correlated Fermi gas 
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We study the BCS-BEC crossover in the strongly correlated regime of an ultra-cold rotating 
two component Fermi gas. Strong correlations are shown to generate an additional long-range 
interaction which results in a modified crossover region compared to the non-rotating situation. 
The two-particle correlation function reveals a smooth crossover between the s-wave paired fermionic 
fractional quantum Hall state and the bosonic Laughlin state. 
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In recent years techniques based on Feshbach scatter- 
ing resonances [l| have enabled the study of pair conden- 
sation in ultra-cold Fermi gases 0, 13, 0| ■ For condensa- 
tion to occur, one can distinguish two distinct physical 
mechanisms: (1) formation of bound pairs of fermionic 
atoms (molecules) which are composite bosons and hence 
undergo Bosc-Einstcin condensation (BEG), and (2) con- 
densation of Bardeen-Goopcr-Schrieffer (BGS) pairs in 
analogy with low temperature superconductivity. In sep- 
arate publications both Eagles and Leggett argued 
that these scenarios were limiting cases of a more general 
theory, the so-called BGS-BEG crossover. 

It was only recently that this crossover phenomenon 
was observed in rotating trap experiments. A vortex lat- 
tice generated in the molecular BEG phase was observed 
to persist into the BGS paired phase as the interaction 
is adiabatically tuned from repulsive to attractive across 
the Feshbach resonance Q- Other experiments, where 
the rotation frequency was increased such that the de- 
generate gas enters the 2D regime, have led to the di- 
rect image of Tkachenko waves Q . In this fast rotation 
regime the effects of strong correlations, such as those re- 
sponsible for the fractional quantum Hall (FQH) effect, 
remain unobserved. It has been predicted that ultra- 
cold atomic systems can be brought into the FQH regime 
by rotating the trap at frequencies close to the trap- 
ping frequency lj 0, In addition, non-rotating 2D 
gases have enabled the observation of the Berezinskii- 
Kosterlitz-Thouless crossover in a trap and in an 

optical lattice [l^ . 

In the present letter we study the BGS-BEG crossover 
in the strongly correlated regime. Usual mean field the- 
ory breaks down as the number of fcrmions becomes com- 
parable to the number of vortices, which is essential to 
reach the FQH regime. We exploit the analogy with a 
FQH system and utilize the Ghern-Simons gauge trans- 
formation which enables us to study superfluidity 
in the strongly correlated regime by considering pairing 
of correlation-free composite fermions. These composite 
particles consist of a fermion and an even number of flux 
quanta. The pairing of fermions in the strongly corre- 



lated regime is then studied by standard BGS mean field 
theory applied to these composite fermions. We show 
that the BGS-BEG crossover in the strongly correlated 
regime can be considered as a crossover between two FQH 
states. 

We consider a two component Fermi system consisting 
of a balanced mixture of fermionic atoms in different hy- 
perfine states represented by ||) and \ [) confined by a 2D 
rotating harmonic trap. In the rotating frame, the Hamil- 
tonian for this system in the FQH regime {Vl — lo ^ Q~) 
in the absence of interactions is given by 
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E 



dr0t(r)^[p + A(r)]'0,(r), (1) 
2 m 



with m the mass of a fermion and A = muyk — mujxy 
is analogous to the vector potential associated with the 
external magnetic field in the electronic FQH effect. The 
operator (/>cr(r) annihilates fermionic atoms with spin a 
at position r. In order to simplify the above Hamiltonian 
we choose to work in a frame where the vector potential 
A is gauged out. This is done by performing the Ghcrn- 
Simons transformation by attaching gauge field Ca{v) = 
-V^Z^cr/ rf^r'eQ/3Pcr(r)(r - - r'p to each bare 

particle resulting in 
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yi^r^^W^ [P + A(i 



)+c(r)f ^.(r), (2) 



where Lp„ is the annihilation operator and pa is the 
density of composite fcrmions of spin tr, and u is the 
filling fraction, which is the ratio of the number of 
atoms to rotational flux quanta. The transformation is 
such that the average gauge field c cancels the exter- 
nal field, thus A(r) + c(r) ~ A(r) -I- c(r) -I- 5c{r) = 
i5c(r). However, we are left with gauge field fiuctua- 
tions which are caused by density fluctuations. If we 
write the interaction part of the above Hamiltonian as 
Hint = (l/2m) Y.. J <Pv^l{v) [2p<5c + ^c^] ^,{v), we see 
that even in the absence of interactions between the bare 
particles, the Ghern-Simons transformation gives rise to 
two- and three-body interactions [l4 1 . 
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The two-body part has been attributed to have impor- 
tant consequences for the formation of pairs in the elec- 
tronic FQH effect • We will neglect the induced three- 
body interaction and consider the induced two-body in- 
teraction y™'^, which has the form 



(3) 



for r < A and which we approximate by V™'^ = for 
r > A. Here A is the typical length scale associated with 
density fluctuations, and we assume that the induced in- 
teraction is washed-out beyond this range. The long- 
range interaction comes on top of the short-range two- 
body interaction, which is dominantly s-wave interaction 
for fcrmions in different spin states. We characterize the 
strength of the interaction via the scattering length a, 
that can be varied by using a Feshbach resonance. This 
resonant short-range interaction remains the same in the 
gauge transformed composite particle picture [l^ . How- 
ever, the additional repulsive long-range interaction can 
strongly modify the resonance properties, and will have 
the effect of lifting up the bound states in the potential 
and change the width of the resonance, resulting in a 
shift and modification of the crossover region. This can 
be treated systematically within the Chcrn-Simons com- 
posite particle picture, where we first need to understand 
the details of this composite particle interaction poten- 
tial. 

The composite fermions experience both a resonant 
short-range and the logarithmic long-range interaction. 
We solve a 2D scattering problem with long-range po- 
tential Eq. ([3]) , where we note that even though the FQH 
effect exists in 2D systems, ultra-cold atomic systems un- 
der extreme rotations are in fact quasi-2D. Quasi here 
means that the confinement in the third dimension is 
strong compared to the remaining two. Hence, the inter- 
action at short range is 3D in nature, and we use the rela- 
tionship between the 2D and 3D scattering length [l^l to 
set a boundary condition at r = 0. This is done by mak- 
ing use of a 2D contact potential [l7| . We solve the 2D 
scattering equations as a function of relative wavenum- 
ber k and the 3D scattering length (osu). For every value 
of asD we can define an energy-dependent 2D scattering 
length a{k), which is related to the scattering phase shift 
6{k) via 



cot S(k) = — (7 
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ka{k) ^ 



(4) 



with 7 the Euler constant. We now characterize the two- 
body interaction strength, which is a result of both in- 
teractions, via a coupling parameter related to the scat- 
tering length calculated for two particles at the Fermi 
energy: 



g2D{kF) = (In 



m kpaikp) 



(5) 
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FIG. 1: (Color online) The 2D interaction strength parameter 
g2D is plotted for different values of A. For this (quasi-2D) 
system a confinement-induced resonance is observed for all 
values of A. For A = .55/o the atomic bound state is "pushed" 
out (for l/{kFa)zD ~ —1.9) by the rotation induced inter- 
action, which inevitably leads to a (narrow) resonance. The 
inset shows the shift and narrowing of the unitarity regime for 
increasing A. The ratio of axial to radial trapping frequency 
is here set to be 10 : 1. 



In Fig.([T]) this coupling parameter is given as function 
of the 3D scattering length, for three different values of 
A. The zero-point oscillations in the third dimension can 
cause density fluctuations which result in gauge field fluc- 
tuations. We believe this will give a natural limit to 
what A may be. On the BCS side we see that beyond 
l/{kpa)^u = —1-9 (for A = .55/o) no two-body bound 
state exists. Since this irnplies that no many-body pair- 
ing instability will occur [l8[ , BCS pairs will be broken. 

In the following, we use a (unitary) Chern-Simons 
transformation to create dressed composite fermions [l9l |. 
This transformation provides a direct relation between 
the wavefunctions of the fermions f) and composite 
fcrmions (^E'cf), 



i<j 



i<3 



X exp 



^|z.|V4-El^^-lV4 



(6) 



where z and ^ are scaled in units of harmonic oscillator 
length lo and represent the complex coordinate of the 
spin up and spin down components respectively. The 
wavefunction of the dressed composite fcrmions is not 
blurred by two-body correlations [20[ . 

We will consider a system at filling fraction = 1/2 
(since this will be the experimentally most accessible 
regime), which means we have transformed fermionic 
atoms to free interacting composite fermions which will 
fill a Fermi sea. Comparison of Eq. ^ to well-known 
paired FQH states, e.g. Haldanc-Rczayi and Moore- 
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Read states, shows that the pairing part of these FQH 
states corresponds to the wavefunction of the compos- 
ite fermions. Now in the presence of some weakly at- 
tractive interaction, for instance caused by an attrac- 
tive atomic interaction, these composites can form BCS- 
like pairs. Hence, we write the Hamiltonian for the 
composite fermion system in the standard BCS form 
and apply (standard) mean field theory. We start 
with the Bogoliubov Hamiltonian in diagonalized form 
HcF = J2k,a ^kll^lka + const.,[2l|, whcrc 7^ ,^ is the 
creation operator for an quasi-particle with energy i?k = 
\/ (ck — m)^ + A^, and ek is the sing le particle kinetic 
energy. Following Randeira et al. [l8| we solve the num- 
ber and gap equation sclf-consistently and find the the 
chemical potential equals /i = 6^?— |£'h|/2, where ep is the 
Fermi energy and Ei, equals the dimer binding energy in 
free space inferred from the energy dependent scattering 
length, both at zero temperature. In the low energy limit 
the gap function becomes a constant A = yj2eF\Eb\. 

Since we consider s-wave spin singlet paired fermions, 
the configuration space first quantized wavefunction for 
2N = N-^ + Ni composite fermions can be written 
as ~ ■!2/(Sii'S22' ...SiVAT') where the anti- 

symmetrization is separately performed over up and 
down spins (the primed and unprimed indexes) [2^ . 



k 



= 'k-(x3-x. 



and ^ i(l±(ek 



/i)/i?k)- Now we are able to determine the form of the 
pairing wavefunctions of the (composite) fermions in the 
BCS and BEC regime. 

To obtain a quantitative view of the crossover be- 
tween the FQH states we consider a system consisting 
of four fermions. On the BCS side we use 'pair' coordi- 
nates {Rem, R, t) where Rem denotes the center of mass 
position, R the distance between the centers of mass 
of the two pairs, and Vjj' = Xj — x^v the interparticlc 
separation of the two fermions forming a pair. These 
four particles have three (equivalent) ways to form pairs 
which all are exphcitly reproduced by Eq. ([6]). Since the 
composite fermions experience a weakly attractive inter- 
action we know \Ei,\ ep and the pair wavefunction 
of the composite fermions for kpr ^ 1 is found to be 
Sjj/ cx sin{kprjj, - n/i)/ y/kpr^j, . 

Via a Feshbach resonance, we tunc the s-wave atomic 
interaction such that the composite fermions experi- 
ence a weakly repulsive interaction, i.e. \Eb\ S> ef. 
The pair wavefunction now describes bosonic molecules 
since S^j/ = {A/2Trep)Ko{Krjji) exactly equals a (deeply 
bound) dimer state in 2D, where k = ik. Advancing to- 
wards the BEC side means the size of the molecule (r) has 
become small as compared to the average inter-particle 
spacing (^ kp^) resulting in 



*F = Sii,S22'(ripr22')'^^^exp ( - R^/A - R, 



cm/ 



(7) 




FIG. 2: The two particle correlation functions G||(r) and 
G|x(r) are shown for N-j — — 100 and ri — Iq. Notice that 
for large r/lo, Gn(r) - G|j.(r) 0. 



consisting of composite bosons, which, form a BEC in the 
presence of the weak repulsive interaction. Transforming 
to 'molecular' coordinates Z, which denote the positions 
of the molecules in complex coordinates, results in 



(8) 



which is a bosonic Laughlin v = 1/% state [23| . The ex- 
ponential contains |Zp/2 terms since the molecular hd.!- 
monic oscillator length equals Iq / \f2. Now we have seen 
a crossover of a paired FQH state a.t i' = 1/2 on the BCS 
side to a bosonic Laughlin i/ — 1/8 state on the BEC 
side for four particles. Note that on the BEC side there 
are half the number of particles and the mass of each 
"elementary" particle (molecule) has been doubled. This 
explains the change of the filling fraction, i.e. the ratio 
of number of particles to the number of flux quanta (in 
this case h/m) from i^=l/2tol/8. 

This crossover between two strongly correlated states 
is believed to be valid for a larger number of fermions. 
To verify this, we have calculated the two-particle corre- 
lation function G(ri — r2) = /■•/ d^r3..rf^rAr|\l/Fp for 
N-^ — Ni = 100 particles using a metropolis Monte- 
Carlo algorithm. Since we know the form of the paired 
FQH state as a function of the BCS coherence length 
(77 = fi^kp /TrmA) we parameterize the Monte-Carlo cal- 
culation by rj. In Fig. [51 we plot both G||(r) and G||(r) 
for tj/Iq = 1. We see that G'||(r) shows a peaked behav- 
ior for small r that is absent in G||(r). At the same time 
for large r, G||(r) — G|x(r) implying the existence 
of a sum rule special to the Haldane Rezayi like state 
of Eq. valid throughout the region of our current 
interest. 

Since in the fc ^ limit, the s-wave T matrix is a 



Being a degree of freedom, the coordinates rii',r22' can 
be effectively integrated out. We have obtained a state 



smooth function of the a^D s-wave scattering length 24| . 
the functional form of the T matrix and hence the gap 
A near the Feshbach resonance will remain unchanged 
hinting a smooth crossover. The crossover behavior of the 
correlation function is clear from Fig.[3l which shows that 
as 7] becomes small compared to Iq, G||(r) gets modified 
continuously and tends towards a limiting form. However 
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FIG. 3; In the strong pairing limit the correlation function 
Gff(r) has a limiting form which equals the G(r) form of the 
1/8 FQH wavefunction Eq.®. 



the most important point to note is that the limiting form 
of Gn(r) is exactly that of the G(r) for the (l/8)-FQH 
state given by the Laughlin form in Eq. ([5]). 

In conclusion, we have shown that the strong corre- 
lations associated with rapid rotations can cause strong 
modifications to the crossover, altering the width and 
position of the crossover. Additionally we have seen the 
fiuctuation induced interaction can be sufficiently strong 
to break BCS pairs and consequently lose superfluidity. 
Using s-wave paired FQH wavefunctions, we have shown 
that the crossover is smooth and the paired FQH state of 
fermions smoothly goes over to 1 /8 bosonic FQH state of 
molecules when one goes across the Feshbach resonance 
such that 7/ <C h- 

A detailed calculation of the crossover physics of this 
region will require the exact nature of the rotation- 
induced long-range interaction Eq. ([3]). Within such 
an treatment for instance Nozieres-Schmitt-Rink calcula- 
tions of the crossover region [2^ can be carried out. Also 
these calculations can be extended to situations with p- 
and d- wave pairing schemes in ultra-cold Fermi gases. 
These scenarios, while having close resemblance with for 
example the 5/2 FQH effect, will be extremely useful and 
will be dealt with in a future publication. 

It is experimentally difficult to reach the FQH regime 
[2^ . A promising techni que is the combination of rota- 
tion with optical lattices [27|, which is for instance able 
to reduce the number of particles per vortex. 

At the same time paired FQH states such as 5/2 are 
known to possess exotic non-abelian quasi-particles exci- 
tations. While existence of non-abelian statistics is the 
basis for topological scheme of implementing quantum 
logic in a quantum computer, the 5/2 state is proved to 
be computationally non-universal. However, there have 
been proposals [28| in which this symptom can be reme- 
died by dynamically tuning-in additional non-topological 
interactions. Dynamic control, while hard in the solid 
state configurations of the FQH effect, transitions be- 
tween different FQH states like the one discussed here 
may be extremely useful for implementing such topolog- 
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